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ABSTRACT 


This  thesis  analyzes  the  performance  of  a  ramjet  with  an 
Oswatisch  inlet  using  a  blunt  centerbody  and  compares  perfor¬ 
mance  to  a  baseline  ramjet  using  an  inlet  with  a  conical 
spike  at  Mach  3.0. 

Inlet  performance  as  a  ratio  of  inlet  lip  to  nose  center- 
body  ratio,  r,/r  ,  is  developed.  The  capture  streamline  for 
each  ratio  is  determined  and  the  coefficient  of  additive  drag 
is  calculated  as  a  function  of  rL/rn.  Setting  thrust  coef¬ 
ficient  equal  to  coefficient  of  drag,  the  performance  of  two 
reunjets  is  determined.  One  ramjet  is  the  baseline  with  a 
spike  inlet;  the  other  ramjet  uses  the  blunt  centerbody. 

Ramjets  and  inlets  are  compared  on  the  basis  of  specific 
fuel  consumption,  excess  thrust  coefficient  and  specific 
thrust.  For  the  ramjet  with  blunt  centerbody,  performance 
parameters  were  calculated  as  a  function  of  inlet  lip  radius 
to  nose  centerbody  radius.  Also  compared  is  the  effect  of 
the  ratio,  r£/rn,  on  relative  detection  range.  For  both 
types  of  ramjets,  the  detection  range  is  reduced  by  approxi¬ 
mately  661.  Performance  of  the  ramjet  with  blunt  nosed 
centerbody  is  severly  handicapped  due  to  high  additive 
drag  and  poor  pressure  recovery.  Specific  fuel  consumption 
is  approximately  50%  greater  for  the  ramjet  with  the  blunt 
centerbody  compared  to  the  ramjet  with  the  spike  inlet. 


TABLE  OF  CONTENTS 


I.  INTRODUCTION . - .  15 

A.  BACKGROUND .  15 

B.  PURPOSE .  19 

II.  FORMULATION  OF  INLET  PERFORMANCE  EQUATIONS 

AND  INLET  FLOW  FIELD .  21 

A.  INLET  DESIGN  OPTIONS . - . r  21 

B.  CALCULATION  OF  STAGNATION  PRESSURE  RATIO 

OF  THE  INLET .  23 

C.  FLOW  FIELD  DESCRIPTION  IN  THE  VICINITY 

OF  A  BLUNT  NOSE .  32 

III.  CALCULATION  OF  RAMJET  PERFORMANCE .  51 

A.  PERFORMANCE  EQUATIONS .  51 

IV.  DISCUSSION  OF  RESULTS .  56 

A.  PERFORMANCE  CRITERION .  56 

B.  RESULTS  OF  CALCULATIONS .  58 

V.  CONCLUSIONS .  69 

APPENDIX  A  -  NASA  Ames  COMPUTER  PROGRAM  DESCRIPTION---  70 

APPENDIX  B  -  MASS  CONTINUITY  METHOD  OF  DETERMINING 

STREAMLINE  SHAPE .  74 

APPENDIX  C  -  ADDITIVE  DRAG  COMPUTATION .  80 

LIST  OF  REFERENCES .  84 

INITIAL  DISTRIBUTION  LIST .  85 


\ 


k. 


5 


LIST  OF  FIGURES 


I- 1  U.S.  Navy  5"/54  Semi  Active  Laser  Guided 

Projectile  (SALGP) .  17 

1 1  - 1  Supersonic  Inlet  With  Blunt  Centerbody -  24 

1 1  -  2  Supersonic  Inlet  With  Conical  Spike -  24 

II - 3  Supersonic  Inlet  (variation) - -  25 

II-4  Supersonic  Inlet  With  Isentropic  Spike .  25 

II -5  Cylindrical  Analog  Of  Busemann  Biplane  Inlet-  26 

II -6  Flow  Regions  Around  A  Blunt  Nosed  Body -  27 

II- 7  Magnitude  Of  Pressure  Losses -  29 

1 1 -  8  Identification  Of  Locations  For  Description 
Of  Pressure  Recovery  Losses  In  A 
Conical  Inlet . 30 

1 1 - 9  Flow  Geometry  Around  The  Blunt  Nose  Showing 

Bow  Shock,  Streamtube  And  Body  Locations -  34 

II- 10  Magnified  View  Of  Geometry  Between  Points 

1  And  2  Of  Figure  II - 9 .  38 

II -11  Geometry  To  Determine  Component  Of  Velocity 
Normal  To  The  Area  For  Integration  Of 
Mass  Flux . . 39 

11-12  Flow  Field  Geometry  For  The  Angular  Method 

Of  Streamtube  Determination . . 44 

11-13  Geometry  For  Additive  Drag  Calculation -  46 

11-14  Geometry  Of  Streamtubes  At  The  Inlet 
Annulus  Used  To  Calculate  Average 
Stagnation  Pressure,  ir^' -  50 

IV- 1  Combustor  Exit  Temperature  As  A  Function 

Of  rL/rn  At  C£  -  CD  »  0.349 .  59 

IV-2  Excess  Thrust  Coefficient  As  A  Function 

Of  r./r . . - . .  60 

Xj  n 

I-V- 3  Specific  Fuel  Consumption  As  A  Function 

Of  r^/rn  For  Ramjet  With  Blunt  Centerbody----  61 


6 


LIST  OF  TABLES 


I I - 1  Capture  Radius,  r  ,  As  A  Function  Of  Lip 

Radius,  r.,  For  r^  Equal  To  1.0  Inch 
At  Body  Station  2  *  1.00 . . 

I I -2  Summary  Of  Radius  Of  Capture  Streamtube, 
Additive  Drag  Coefficient,  And  Mass  Flow 
Ratio . . — . — 

IV- 1  Comparison  Of  Spike  Inlet  Versus  Two  Blunt 

Nose  Inlets - . - 

A-l  Reference  Values  For  OGIVE  and  IMPLCBO 

Programs - 

B-l  Variables  For  Capture  Streamtube  Computer 

Program - - - 

B-2  Capture  Streamtube  Program  Listing - 

C-l  Variables  For  Coefficient  Of  Additive 

Drag  Computer  Program- . — 

C-2  Coefficient  Of  Additive  Drag  Program 

Listing . 


ll 


8 


LIST  OF  SYMBOLS 


Symbol 

Explanation 

Units 

a 

Speed  of  sound 

ft/sec 

A 

Area  (with  subscript) 

ft2 

Ar 

Reference  area  which  is  maximum 
cross  sectional  area  for  pro¬ 
jectile 

ft2 

Cdad 

Coefficient  of  additive  drag 

CP 

Coefficient  of  specific  heat 
capacity  at  constant  pressure 

BTU/lb  • 
m 

Cf 

Coefficient  of  thrust 

CD 

Coefficient  of  airframe  drag 

C£e 

Coefficient  of  excess  thrust 

D 

a 

Additive  drag 

f 

Fuel  air  ratio  (ra^/mo) 

F 

Thrust 

lbf 

g 

Gravitational  constant 

ft/sec 

h 

Heating  value  of  hydrocarbon 

BTU/lbm 

m 

Mass  flow  rate 

slug/sec 

lb^/sec 

■s 

Mass  flow  rate  in  the  capture 
streamtube  at  the  bow  shock 

lb  /sec 
m 

“5 

Mass  flow  rate  (reference)  equal 
to  mass  flow  rate  in  a  circular 
streamtube  of  radius  r^ 

lb  /sec 
m 

p 

Pressure  (static  without 
subscript) 

lb/ft2 

Q 

Non  dimensional  free  stream 

velocity  defined  as  u  /v 

s  m 

9 


Units 


Symbol 


r 


s 


r 


L 


R 


r 


S 

SF 

SFC 

Sr 

T 

u 


v 


X 

X 


Explanation 

2 

Free  stream  dynamic  pressure  lb^/ft 

Nose  radius  of  inlet  centerbody;  in 
nose  radius  of  seeker 

Radius  of  the  capture  streamline  in 
at  the  intersection  with  the  bow 
shock  wave 

Outer  radius  of  the  inlet  annulus  in 


Distance  normal  to  projectile  in 
axis  for  region  2  coordinate 
system 


Distance  from  projectile  axis 
using  cylindrical  coordinates; 
r  is  identical  to  Y  used  in  the 
computer  printout 

Capture  streamtube  area 

Specific  thrust 

Specific  fuel  consumption 

Quantity  defined  by  equation  44 

Temperature 

Component  of  velocity  in  the  Z 
or  X  direction 


in 


lbf/(lbm/sec) 

(lbm/hr)/lbf 


°R 


ft/sec 


Maximum  velocity  obtained  when 
a  gas  with  stagnation  speed  of 
sound  (a.)  is  expanded  into  a 
vacuum;  v]j*(2/(7-l))1'2  at  ft/sec 

Component  of  velocity  normal  to 
projectile  axis  in  region  1; 
for  ramjet  performance 

Velocity  of  air  at  designated  ft/sec 

station  of  ramjet 

Ratio  as  defined  in  equation  30 

Distance  parallel  to  projectile  in 
axis  for  region  1  coordinate 
system 


10 


Units 


Symbol 

Z 


a 

8 

Y 

n 

9 


Explanation 

Distance  parallel  to  projectile 
axis  for  region  2  coordinate 
system 

Distance  perpendicular  to 
projectile  axis  for  region  1 
coordinate  system 

Angle  of  attack 

Angle  measured  from  reference 
direction  to  the  local  velocity 
flow  vector 

Ratio  of  heat  capacities 
Efficiency 

Angle  of  streamtube  between 
two  data  points 

Ratio  of  downstream  stagnation 
pressure  to  upstream  stagnation 
pressure  (with  subscript);  Ratio 
of  circle  circumference  to  circle 
diameter  (without  subscript) 


Density 


in 

in 

degrees 


degrees 


11 


Subscripts 

a  Additive  drag 

b  Burner 

d  Diffuser 

e  Exit  station  of  ramjet  engine 

f  Force  (with  lb) 

i  Integer  identifying  points  along  capture 

streamline 

j  Integer  identifying  individual  streamtubes  within 

the  capture  streamline 

L  Lip  of  inlet  annulus 

m  Mass  (with  lb) 

n  Normal  to  surface;  nozzle  with  ramjet;  nose  of 

proj  ectile 

o  Inlet  station  of  ramjet  engine 

s  Free  stream  values  at  the  bow  shock  on  the 

upstream  side 

t  Stagnation  values 

1  Designates  values  from  region  1  as  defined  in 
Figure  1 1 - 6 

2  Designates  values  from  region  2  as  defined  in 
Figure  I I -6 

3  Combustor  exit  station 

a  Values  along  a  constant  line  as  indicated  in 

Figure  1 1 - 9 

00  Free  stream  values 


12 


Superscripts 


(  ) •  Values  tabulated  in  NASA  Ames  computer  printout 

(non  dimensional) 

(  ~ )  Average  value 


ACKNOWLEDGMENT 


The  author  wishes  to  express  his  thanks  to  Dr.  Paul 
Kutler  and  Dr.  Dennis  Chaus se  of  NASA  Ames  Computational 
Fluid  Dynamics  Branch  for  their  assistance  in  the  calcula¬ 
tion  of  the  flow  field  over  the  blunt  noses  body.  Their 
assistance  aided  this  work  immeasurably. 

Also,  the  author  wishes  to  express  his  deepest 
gratitude  to  Dr.  Allen  E.  Fuhs  for  his  patience,  guidance 
and  time.  Without  him,  this  thesis  would  not  have  been 
possible. 


14 


I.  INTRODUCTION 


A.  BACKGROUND 

The  rapid  development  and  operational  deployment  of  long 
range  antiship  cruise  missiles  has  generated  demanding 
requirements  for  the  anti-ship  missile  defense  (ASMD)  mis¬ 
sion.  The  threat  has  necessitated  that  the  Navy  develop  the 
concept  of  a  "defense  in  depth"  for  the  defensive  doctrine 
of  a  carrier  battle  group  (CBG) .  The  doctrine  emphasizes 
the  utilization  of  all  available  weapon  systems  in  a  layered 
defense  to  defeat  hostile  targets.  Conceivably,  an  anti- 
ship  missile  (ASM)  could  survive  attacks  from  extremely  long 
range  missiles,  the  E-2C  and  F-14  combat  air  patrol  (CAP) 
team.  Standard  ER  and  Standard  MR  missiles,  conventional  gun 
ordnance,  point  defense  weapons  and  close  in  defense  weapons, 
such  as  Phalanx,  prior  to  reaching  a  high  value  unit.  Since 
enemy  doctrine  is  to  fire  the  ASM  in  such  numbers  as  to 
overwhelm  the  Command  and  Control  capability  of  the  CBG 
and  of  an  individual  vessel,  the  accent  of  new  ASMD  weapons 
systems  is  on  quick  response  time,  short  time  of  flight  and 
a  high  degree  of  accuracy  against  a  manuvering  ASM. 

Although  a  single  saturation  raid  could  possibly  be 
defeated,  what  about  a  second  or  third  raid  in  a  single  day? 
Will  the  CBG  have  enough  missiles  to  supplement  airborne  CAP 
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and  still  destroy  the  ASM  threat  30  to  40  miles  from  the 
high  value  unit? 

One  solution  to  this  vexing  problem  is  to  modify  the 
5"/54  Semi-Active  Laser  Guided  Projectile  (SALGP)  which  is 
illustrated  in  Figure  1-1.  The  SALGP  round  is  currently  un¬ 
dergoing  operational  evaluation  (OPEVAL)  prior  to  fleet  wide 
introduction  [Ref.  1]  and  consists  of  a  five  inch  projectile 
with  a  laser  seeker  on  the  projectile  nose.  The  round  is 
designed  to  "home  in"  on  laser  energy  reflected  from  a  target 
by  a  laser  illuminator.  The  modification  to  the  SALGP  round 
consists  of  replacing  the  laser  seeker  head  with  an  infrared 
(IR)  seeker  and  adding  a  solid  fuel  ramjet  to  power  the  pro¬ 
jectile.  The  ramjet  would  provide  constant  speed  over  a 
longer  portion  of  the  flight  path  and  provide  energy  in  the 
terminal  phase  of  flight  to  allow  for  manuvering  to  hit  a 
manuvering  target.  The  modified  SALGP  is  identified  as  a 
gun  launched  missile  (GLM) .  The  modifications  would  allow 
the  GLM  to  become  a  fire  and  forget  round  requiring  no  ter¬ 
minal  guidance  from  the  ship. 

Although  not  helping  in  the  long  range  defense  of  a  CBG, 
the  GLM  would  complement  the  point  defense  and  close  in 
defense  weapon  systems  to  the  ships  in  the  U.S.  Navy  today. 

The  GLM  will  allow  a  naval  vessel  without  a  missile  sys¬ 
tem  and  with  a  properly  configured  gun  system  to  increase 
its  anti-air  warfare  (AAW)  capability  in  a  cost  effective 
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way.  More  GLM  rounds  can  be  carried  in  existing  gun  maga¬ 
zines  than  missile  rounds  can  be  carried  in  current  missile 
magazines.  Within  a  ten  nautical  mile  zone  around  a  ship, 
the  GLM  will  have  a  faster  response  time,  a  shorter  time  of 
flight  and  be  as  accurate  as  guided  missiles  in  the  current 
Navy  inventory.  Fast  response  and  short  time  of  flight  are 
due  to  the  initial  velocity  imparted  to  the  GLM  by  the  5”/ 54 
gun.  The  shape  of  the  nose  of  the  GLM  round  will  be  con¬ 
strained  by  the  use  of  the  hemispherical  lense  surface  for 
the  infrared  seeker  in  the  projectile  nose.  This  blunt  nose, 
although  optimized  for  the  sensor  optics,  is  not  necessarily 
the  best  aerodynamic  shape. 

Although  the  SALGP  has  the  potential  to  be  highly  effec¬ 
tive  in  the  Naval  Gunfire  Support  (NGFS)  and  in  the  Surface 
to  Surface  Warfare  (SUW)  missions,  it  was  not  designed  ini¬ 
tially  for  the  Air  Warfare  (AAW)  or  Anti-ship  Missile 
Defense  (ASMD)  missions.  Consider  a  GLM  based  on  SALGP 
hardware  to  the  maximum  extent  possible.  To  become  success¬ 
ful  in  the  AAW/ASM D  mission  areas,  it  is  essential  to  im¬ 
prove  the  overall  SALGP  performance  and  capability.  One  area 
to  achieve  easily  measured  performance  gains  is  in  the 
method  of  propelling  the  SALGP. 

One  propulsion  proposal  that  decreases  mission  flight 
time  and  increases  mission  effectiveness  is  to  replace  the 
solid  propellant  rocket  on  the  SALGP  round  with  that  of  a 


ramjet  with  either  liquid  or  solid  fuel.  Obviously,  to  keep 
total  system  developement  and  procurement  costs  low,  and 
therefore  warrant  the  use  of  a  ramjet  in  a  cost  effective 
analysis,  the  SALGP  round  should  be  changed  as  little'  as 
possible  in  adding  the  ramjet.  Effectively,  the  designer  is 
constrained  to  use  a  blunt  nosed  centerbody  to  house  the  IR 
seeker  and  sensor  package.  The  remainder  of  the  projectile 
is  constrained  by  the  dimensions  of  the  handling  equipment 
in  the  MK  45  mount  and  Department  of  the  Navy  Specifications 
[Ref.  2]. 

B .  PURPOSE 

Owe  purpose  of  this  thesis  is  to  investigate  the  effect 
that  a  v.  tunt  nosed  centerbody  has  on  the  pressure  recovery 
of  an  inlet.  The  pressure  recovery  is  an  important  indica- 
:or  of  ramjet  performance,  for  poor  pressure  recovery  will 
reduce  the  flight  envelope  over  which  the  ramjet  operates. 

Although  drag  is  large,  the  blunt  nosed  body  was  aerody- 
namically  feasible  for  a  rocket  guided  projectile,  since  the 
operation  of  the  rocket  engine  does  not  require  air  to  be 
"brought  onboard".  The  use  of  a  ramjet  may  prove  inadvise- 
able  or  infeasible  due  to  aerodynamics  of  the  flow  around 
the  blunt  nose  body  enroute  to  the  inlet  annulus.  The  pres¬ 
sure  recovery,  or  lack  of,  may  also  require  the  body  shape 
surrounding  the  IR  sensor  system  be  reconfigured,  within 
existing  system  constraints,  to  allow  for  ramjet  operation. 
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The  aerodynamic  model  used  was  that  of  a  hemispherical  body 
with  an  attached  cylindrical  body  of  constant  diameter.  The 
flow  was  calculated  for  a  cylinder  length  equal  to  four 
times  the  nose  radius.  Distances  involved  were  normalized 
to  the  nose  radius,  r  .  This  is  the  same  as  the  radius  of 
the  hemisphere. 

The  National  Aeronautics  and  Space  Administration  Ames 
Research  Center  (NASA  Ames)  calculated  the  flow  over  a  hem¬ 
ispherical  body  with  an  attached  cylinder  over  the  range  of 
Mach  numbers  from  1.8  to  3.4,  in  steps  of  0.2,  at  zero  angle 
of  attack  and  at  Mach  3.4  with  an  angle  of  attack  (oi)  equal 
to  10  degrees.  The  calculations  provide  flow  field  data 
which  included  shock  shape,  sonic  line,  velocity  components, 
pressure,  density,  entropy  and  similar  flow  parameters.  The 
calculations  were  done  by  computer  code  contained  in  two 
separate  programs.  For  a  further  discussion  of  the  programs 
and  the  various  inputs  see  Appendix  A. 

Section  II  reports  on  inlet  performance  and  the  inlet 
flow  field.  Section  III  reports  on  ramjet  performance  for  a 
blunt  nosed  centerbody  for  the  inlet  and  compares  the  per¬ 
formance  of  a  ramjet  with  an  isentropic  spike  inlet. 
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II.  FORMULATION  OF  INLET  PERFORMANCE  EQUATIONS 

- AND'  INLETTLC’W  HELP - - 

A.  INLET  DESIGN  OPTIONS 

In  developing  a  quantitative  basis  for  comparison  of  var¬ 
ious  inlet  designs,  several  items  of  information  are  needed. 

First,  information  needed  is  the  drag  at  the  different 
Mach  numbers  and  a  correlation  of  drag  versus  Mach  number. 

The  blunt  nosed  centerbody  produces  a  detached  shock  that 
stands  off  a  distance  from  the  body.  A  conical  inlet  pro¬ 
duces  an  attached  oblique  shock. 

Second,  the  on-design  performance  of  the  inlet  design 
must  be  evaluated.  In  particular,  the  pressure  recovery  of 
the  inlet  is  important.  What  happens  to  this  recovery  ratio 
when  the  position  of  the  annulus  is  moved  axially  along  the 
body?  If  the  location  is  fixed  and  the  radius  of  the  inlet 
lip  varied  what  happens  to  the  inlet  performance?  A  compar¬ 
ison  must  be  made  of  the  relative  ranges  attained  by  a  ram¬ 
jet  with  an  isentropic  spike  and  one  with  a  blunt  nosed 
centerbody  must  be  made.  Also  of  concern  is  how  sensitive 
the  ramjet  performance  is  to  angle  of  attack  variations. 

Third,  off-design  performance  of  the  ramjet  is  of  inter¬ 
est  in  several  aspects.  First,  how  does  the  pressure  recovery 
of  the  inlet  vary  at  different  Mach  numbers?  Ramjet  perfor¬ 
mance  is  sensitive  to  the  inlet  pressure  recovery.  Generally 
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ramjet  performance  tends  to  improve  initially  at  higher  Mach 
numbers.  Second,  how  does  this  affect  projectile  drag  and 
what  effect,  if  any,  is  transmitted  to  the  projectile  range? 
Third,  what  is  the  operational  region  over  which  the  inlet 
and  ramjet  will  continue  to  operate?  The  lower  performance 
limit  of  the  inlet  is  that  Mach  number  where  the  normal  shock 
within  the  inlet  is  expelled.  This  condition  is  refered  to 
as  subcritical  ramjet  operation  and  can  lead  to  a  condition 
called  "buzz".  The  upper  performance  limit  of  the  inlet  is 
the  maximum  flight  speed  at  which  adequate  pressure  recovery 
can  be  attained.  Last,  how  does  the  angle  of  attack  affect 
off  design  performance? 

The  desire  that  the  ramjet  powered  GLM  fit  within  the 
existing  confines  of  the  5"/54  MK  45  gun  limits  the  inlet 
design  to  that  of  the  axisymmetric  variety.  In  theory, 
"pop-up"  inlets  or  scoops  could  be  used.  However,  the  ability 
of  the  inlets  to  survive  a  gun  launch  of  8000  g's  and  func¬ 
tion  reliably  increases  the  complexity  of  the  round;  variable 
geometry  inlets  are  not  a  good  use  of  the  remaining  available 
space  in  the  projectile.  The  small  diameter  of  the  projectile 
coupled  with  the  nose  mounted  laser  seeker,  requires  intel¬ 
ligent  use  of  available  space. 

The  axisymmetric  inlets  which  could  be  utilized  are 
1)  blunt  nosed  centerbody,  2)  conical  nose  and  variations, 

3)  an  unconstrained  design,  an  isentropic  spike  for  example, 
and  4)  a  cylindrical  analog  of  the  Busemann  Biplane. 


The  blunt  nose  shown  in  Figure  I I - 1  is  currently  in  use 
on  the  SALGP  projectile,  without  the  annulus.  As  modeled  by 
the  simulation,  the  inlet  consists  of  a  hemispherical  nose 
connected  to  a  constant  diameter  body  extending  beyond  the 
inlet  lip.  The  conical  nose  in  Figure  1 1  - 2  uses  a  cone  of 
fixed  angle  to  focus  the  shock  wave  at  the  inlet  lip  when 
at  the  design  Mach  number.  A  variation  is  shown  in  Figure 
1 1 - 3 .  Shown  in  Figure  1 1 - 4  is  the  unconstrained  design. 

The  design  would  consist  of  an  isentropic  spike.  Figure 
1 1 - 5  shows  the  cylindrical  analog  of  the  Busemann  Biplane. 

The  inlet  "swallows"  the  oblique  shock  wave  formed  by  the 
nose  at  the  design  Mach  number  and  has  a  body  of  constant 
diameter.  The  advantage  is  the  elimination  of  external  wave 
drag. 

B.  CALCULATION  OF  STAGNATION  PRESSURE  RATIO  OF  THE  INLET 

For  calculation  purposes,  the  flow  around  the  blunt  nosed 
body  is  divided  into  two  segments  as  shown  in  Figure  II-6. 
Region  1  is  the  region  of  flow  from  upstream  infinity  to  the 
shoulder  of  the  centerbody.  The  centerbody  shoulder  is 
located  in  plane  S  in  Figure  II-6.  Region  1  is  characterized 
by  flow  in  both  the  longitudinal  and  radial  directions. 

Region  2  is  the  region  of  flow  downstream  of  plane  S.  The 
majority  of  the  flow  is  in  the  longitudinal  direction  with  a 
relatively  small  radial  component.  The  origin  for  coordi¬ 
nate  systems  for  regions  1  and  2  is  the  same  as  illustrated 
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plane  S 


Figure  1 1 - 6  Flow  Regions  Around  A  Blunt  Nosed  Body 


2 


in  Figure  II -6.  The  coordinates  for  region  1  are  X  and  Y, 
which  are  the  symbols  used  in  the  NASA  Ames  computer  printout. 
The  coordinates  for  region  2  are  R  and  Z,  also  illustrated 
in  Figure  II -6.  Each  of  the  regions  has  different  reference 
values  and  symbols  for  the  flow  properties.  Appendix  A  lists 
the  reference  values  and  symbols. 

In  a  real  ramjet,  the  sizing  of  the  inlet  capture  area, 
throat  area  and  inlet  area  can  dramatically  vary  the  perfor¬ 
mance  at  different  Mach  numbers.  The  importance  of  the  effect 
the  inlet  plays  on  performance  is  paramount. 

Figure  II- 7  shows  the  magnitude  of  pressure  losses  as  a 
fraction  of  total  pressure  recovery.  As  can  be  clearly  seen, 
the  majority  of  loss  is  due  to  boundary  layer  growth.  How¬ 
ever,  as  the  flight  Mach  number  increases,  note  that  all  the 
losses  increase  in  magnitude. 

For  a  conical  inlet,  the  oblique  shock  losses  can  readily 
be  calculated.  The  overall  pressure  recovery  value  of  0.72, 
for  Mach  3.0  flow,  can  be  obtained  from  Figure  1 1 - 7 . 

For  the  blunt  nosed  body,  pt/ptoo  is  tabulated  in  the  NASA 
Ames  printout.  This  corresponds  to  P^/Ptl*  oblique 
shock  losses,  in  Figure  1 1 - 8 .  To  account  for  the  varying 
pressure  recovery  values  within  the  annulus,  a  mass  weighted 
average  will  be  used.  This  concept  will  be  discussed  in 
Section  II-C-3. 
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Figure  II- 8  Identification  Of  Locations  For 
Description  Of  Pressure  Recovery 
Losses  In  A  Conical  Inlet 
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Total  pressure  ratio,  or  pressure  recovery,  is  defined 


as  From  Figure  1 1  - 8 , 


£ts  hi  hi 
Ptl  *  Ptl  Pt2 


(1) 


By  defining  as  equal  to  P^^Ptl*  t^ie  Pressure  recovery 
becomes 


*d 


t2) 


The  value  of  Pt$/Pt2  can  determined  from  Figure  1 1 - 7 . 
The  value  of  Pt5/Pt2  represents  the  accumulation  of  losses 
in  an  inlet  due  to  the  subsonic  diffuser  loss,  p^/p^; 
transonic  shock  and  viscous  losses,  Pt4/Pt3»  an<i  pressure 
loss  due  to  boundary  layer  growth,  Ptj/Pt2*  The  value  of 
Pt5/pt2  *s  obtained  from  Figure  1 1  -  7  as  a  ratio  of  overall 
inlet  pressure  recovery  to  oblique  shock  losses.  For  Mach 
3 . 0  flow 


hi 

pt2 


0.72 

TTT97 


0.7423 


The  value  of  was  calculated  using  the  flow  field  calcu¬ 
lations  in  the  NASA  Ames  printout.  The  overall  effect  of  the 
reduced  recovery  ratio  is  that  the  engine  efficiency  and  the 
pressure  ratio  at  the  nozzle  exit  is  changed.  In  turn,  the 


change  in  the  pressure  ratio  may  cause  over  or  under  expanded 
flow  from  the  exit  nozzle  [Ref.  3], 

C.  FLOW  FIELD  DESCRIPTION  IN  THE  VICINITY  OF  A  BLUNT  NOSE 
For  the  blunt  nosed  centerbody  at  zero  angle  of  attack, 
the  flow  field  is  axisymmetric  with  respect  to  the  centerline 
of  the  projectile.  The  calculation  of  pressure  recovery  and 
ramjet  performance  is  greatly  simplified. 

1 .  Determination  Of  Capture  Streamline 

Inherent  to  the  analysis  of  ramjet  performance  is  the 
determination  of  the  mass  flow  rate  of  air  at  the  annulus  of 
the  inlet,  m^.  After  determining  the  mass  flow  rate,  mL, 
the  streamline  coordinates  will  be  determined  at  various 
points  by  equating  the  local  mass  flow  to  the  mass  flow  at 
the  lip.  From  these  coordinates,  the  additive  drag  and 
coefficient  of  additive  drag  for  the  projectile  at  a  parti¬ 
cular  Mach  number  will  be  determined.  Two  approaches  for 
calculating  the  streamtube  points  were  tried  and  will  be 
discussed  here. 

a.  Mass  Flow  Method 

The  first  method  involves  mass  flow  analysis. 

The  mass  flow  through  any  cross  section  of  the  particular 
streamtube  will  equal  that  of  the  mass  flow  through  the 
annulus.  Of  particular  interest  is  the  radius  of  the  capture 
area,  rg,  on  the  bow  shock  wave.  The  air  captured  by  this 
area  goes  into  the  inlet.  The  remainder  of  the  air  spills 
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past  the  inlet  and  affects  the  performance  of  the  projectile 
by  adding  to  the  total  vehicle  drag  with  a  term  called  addi¬ 


tive  drag,  D  . 

a 

The  mass  flow  at  the  lip  is  given  as 


m 


L 


dA 


(3) 


Where  p^  is  the  density,  dA  the  elemental  area  and  the 
velocity  component  normal  to  that  area.  From  the  continuity 
equation,  the  mass  flow  at  the  shock  front  is 


P^r 
s  s 


u. 


(4) 


Refer  to  Figure  1 1 - 9 .  By  equating  equations  3  and  4 


p  ffr  u 
s  s  s 


2nrL  drL 


CS) 


By  normalizing  the  distances  involved  to  the  nose  radius  and 
by  dividing  by  pM  and  v  ,  to  allow  use  of  the  NASA  Ames 
printout  values  in  region  2  ,  the  equations  become 


7  Pt°°  /’P  '  '  j  ' 

2  PT  J°2  U2  r  dr 


body 


(6) 


Plane  S 


Figure  I I -9  Flow  Geometry  Around  The  Blunt  Nose 
Showing  Bow  Shock,  Streamtube  And 
Body  Locations 
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where  u. 


‘2  ul/vm:  H2  "L' "t®’  -  ’  '•  n 

values  are  the  tabulated  values  from  the  NASA  Ames  printout 

and  the  subscript  refers  to  the  applicable  flow  region  de¬ 
fined  in  Figure  I I -6.  Q 
becomes 


-  pL/pt0O;  and  r 


r/r  .  The  primed 


is  defined  as  u  /v  ,  and  equation  6 
s  m 


r2 

s 

T 

rn 


2  p 


•Hp 
1 2  “2 
'body 


^  /* Po’  u -i  r  dr 

00  A, 


(7) 


The  integrand  in  equation  7  was  plotted  and  the  area  under 
the  curve  calculated  for  the  desired  lip  radius.  Since 
r  =  1.0  inches,  Q  =  0.80174,  p  /p..  =  0.07623  at  Mach  3  and 

the  integral  from  body  to  lip  is  0.00805,  rg  can  be  determined. 
For  Mach  3  airflow,  with  r^  *  1.2 

rs  =  ( T578UT7  <f)T0  •  GTdl  3)  CO. 00805) 


Hence,  r  *  0.5133. 

*  s 

Table  II-l  shows  rg  for  Mach  3.0  and  a  nose  centerbody  radii 
equal  to  1.00  inches.  The  maximum  value  of  the  lip  radius 
is  1.854  inches.  This  is  the  size  imposed  by  the  constraints 
of  the  5"/54  handling  system.  Now  the  streamline  locations 
at  different  flow  field  points  can  be  calculated.  The  compu¬ 
ter  program  used  is  described  in  Appendix  B.  The  flow  geometry 
is  shown  in  Figure  I I -9. 


Table  I I - 1 

Capture  Radius,  rg,  As  A  Function  Of  Lip  Radius,  ,  For  r^ 
Equal  To  1.0  Inch  At  Station  Z=*1.00 


rL  (inches) 


1.1000 

1.2000 

1.3000 

1.4000 


1.5000 


1.6000 

1.7000 

1.8000 

1.8540 


rg  (inches) 


0.3821 

0.5883 

0.7765 

0.9610 


1.1474 


1.3353 

1.5231 

1.7066 

1.8086 


The  mass  flow  rate  across  a  line  of  constant 


angle,  a,  is  given  by  the  following  equation 


m 


a 


dS 


(8) 


where  b  refers  to  the  body  and  cs  refers  to  the  capture  stream¬ 
line  in  Figure  1 1 -  9 .  From  Figure  11-10,  an  element  of  area 
on  the  line  of  constant  a  is  given  as 


dS  =  2irrdl 


(9) 


and  dl  can  be  related  to  the  angle  a  by 


dl 


dr 

sina 


(10) 


Hence,  dS  in  terms  of  dr  becomes 


dS 


2irr 


dr 

sina 


(ID 


From  Figure  II-ll,  an  angle  3  is  defined  so  that 


3 


tan  * (v/u) 


(12) 
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Figure  I I -10  Magnified  View  Of  Geometry  Between 
Points  1  And  2  Of  Figure  1 1 - 9 
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Figure  11-11  Geometry  To  Determine  Component  Of  Velocity 
Normal  To  The  Area  For  Integration  Of  Mass 
Flux 


Hence,  the  normal  component  of  the  velocity  becomes 


=  /  u2  +  v2  cos  C  tt/ 2  -  ex  -  3)  (13) 

With  the  substitution  of  a  trigonometric  identity  for  the 
cosine  term,  equation  13  can  be  simplified  to 

- 2 - 

u^  =  /  u  +  v  sin(a+8)  (14) 

The  geometry  of  Figure  I I- 11  gives 


cos  8 


and 


sin  3 


(15) 


(16) 


Using  equations  15  and  16  and  substituting  into  equation  14 
leaves 
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In  front  of  the  projectile  shoulder,  equation  21  becomes 


2 


2 


III  I 

r  (u^  +  cotoOdr 


(22) 


For  region  2  ,  which  is  after  the  shoulder,  due  to  the  dif¬ 
ferent  reference  values  equation  20  becomes 

r  t  2p  C  '  '  •  *  ' 

C—)^  *  Q— 5-  /  P2  r  (u2  +  v  cota) dr  (23) 

b 

where  in  equations  22  and  23,  the  primed  values  are  taken 
from  the  NASA  Ames  printout.  In  region  2  ,  a  equals  90°  and 
cotangent  of  a  is  zero. 

b.  Angular  Method 

The  second  method  of  streamtube  calculation  is 
by  the  angular  method.  A  linear  variation  of  the  known  flow 
field  properties  over  small  distances  is  assumed.  The  method 
starts  at  the  lip  and  progressively  moves  forward  along  the 
body  determining  streamtube  position. 

Initially,  the  local  velocity  vector  at  the 
various  flow  field  points  must  be  plotted  in  polar  notation. 

By  averaging  the  velocity  flow  angles  at  two  adjacent  stations, 
an  average  angular  value  can  be  determined.  As  an  example, 
for  a  lip  radius  of  1.2 
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where  03  ■  local  flow  field  velocity  angle  at  station  3 

@4  =  local  flow  field  velocity  angle  at  station  4 

0,,=  average  angle  between  station  3  and  4  stream 
tube  locations 

04,  the  flow  field  angle  at  the  lip,  is  determined 
from  linear  interpolation  of  the  known  flow  field  angles.  By 
varying  043,  with  0^  known,  the  value  of  03  is  determined. 

As  a  function  of  843. 

04  =  2043  -  03  (25) 

043  is  then  plotted  as  a  function  of  the  radius  of  the 
streamtube,  R3  in  this  case. 

By  studying  Figure  11-12,  the  geometry  of  the 
field  also  determined  043- 

R  -  R 

tan  943  ’  [  r?—  -ZT  1  w 

4  3 

In  the  above  equation,  all  quantities  with  the  exception  of 
R3  are  known.  043  is  then  plotted  versus  R3. 


The  intersection  of  these  two  curves  determines 
Rj  and  0^,  and  hence  0^  can  be  determined  by  interpolation. 
In  this  manner,  the  streamtube  points  are  calculated  for  the 
desired  lip  radius.  This  method  produces  comparable  results 
to  the  mass  flow  analysis  but  is  not  accurate  enough  pro¬ 
ceeding  past  the  nose  of  the  projectile.  Additionally,  the 
process  is  extremely  slow  and  complicated. 

2 .  Determination  Of  Additive  Drag 

The  additive  drag  is  calculated  b>  raputer  program 

discussed  in  Appendix  C.  Assuming  that  the  points  defining 
the  edge  of  the  streamtube  determined  by  mass  flow  analysis 
can  be  connected  by  straight  line  segments,  the  equation  for 
additive  drag  can  be  written  as 


Da  *  EfPi  •  PJai  sin9i 


(27) 


where  p^  equals  the  average  pressure  at  the  data  points  i 
and  i+1;  0^  is  the  angle  of  the  traight  line  segment  between 
the  points  i  and  i+1  measured  from  the  projectile  axis;  and 
A^  is  the  surface  area  of  the  cone  between  the  i  and  i+1  data 
points.  The  geometry  is  shown  in  Figure  11-13. 

The  coefficient  of  additive  drag  (C^ad^  can  be  de¬ 
termined  by  dividing  equation  27  by  q^  and  reference  area,  Ar, 
2  2  P.4±  -  1)  At  sinet 

*  00 


Cdad 


YP. 


(28) 
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Figure  11-13  Geometry  For  Additive  Drag  Calculation 
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•  t  -  <i 


The  reference  area,  Ar,  is  taken  to  be  the  base  of  the  five 
inch  projectile.  Simplification  of  equation  (28)  leads  to 


Jdad 


ymq 


tpi  -  i) 

ST  qo 


Ai 

JT  sin8i 


(29) 


Equation  (29)  is  the  equation  programmed  in  Appendix  C. 

The  program  utilizes  the  points  generated  by  the  pro¬ 
gram  in  Appendix  B  to  calculate  C^a(j.  However,  it  is  impor¬ 
tant  to  scale  the  nose  capture  radius  from  a  centerbody  radius 
of  one  inch  to  the  actual  centerbody  radius  on  the  projectile. 
Previously  r^  was  fixed  at  1.854  inches.  A  ratio  between 
the  centerbody  and  lip  radii  is  defined  as  x  *  r^/r^. 

The  lip  radius,  rL,  is  f:xed  by  the  dimensions  of  the 
handling  equipment,  and  rn  can  be  determined  as  a  function  of 
the  ratio  x.  Table  1 1  -  2  displays  the  values  of  x,  r  , 

and  mass  flow  ratio.  Mass  flow  ratio  is  defined  as  m  /m°. 

s  s 

The  symbols  ms  and  ra°  are  defined  in  Table  1 1 - 2 . 

3 .  Average  Stagnation  Pressure  Ratio  At  The  Inlet  Lip 
The  stagnation  pressure  at  the  annulus  of  the  inlet 
is  a  function  of  radius.  A  mass  weighted  average  stagnation 
pressure  was  calculated.  The  average  pressure  recovery  will 
be  multiplied  by  a  factor  of  Pts/Pt2  t0  obtain  total  pressure 
recovery  as  discussed  in  section  II. 
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Table  1 1 -2 


Summary  Of  Radius  Of  Capture  Streamtube,  Additive  Drag  Coef 

ficient.  And  Mass  Flow  Ratio 


1.1 

0.3821 

1.6855 

0.6440 

0.3726 

0.121 

1.2 

0.5883 

1.5450 

0.9089 

0.2956 

0.240 

1.3 

0.7765 

1.4262 

1.1074 

0.2299 

0.357 

1.4 

0.9610 

1.3243 

1.2727 

0.1775 

0.471 

1.5 

1.1474 

1.2360 

1.4181 

0.1319 

0.595 

1.6 

1.3353 

1.1588 

1.5473 

0.0928 

0.697 

1.7 

1.5231 

1.0906 

1.6611 

0.0588 

0.803 

1.8 

1.7066 

1.0300 

1.7578 

0.0304 

0.899 

1.854 

1.8086 

1.0000 

1.8086 

0.0179 

0.952 

rn  is 
s 

defined  by 

equation  (4)  ; 

m°  is  obtained  from 
s 

equation 

(4)  with  rg  replaced  by  rL 


The  mass  weighted  average  of  pressure  recovery  is 


given  by 


i 


x 

I 

j=l 


Pt  m. 

C— )  • 

Pt»  J  mL 


(30) 


where  the  subscript  j  refers  to  the  jth  streamtube  within 
the  inlet,  m.  is  the  mass  flow  of  the  jth  streamtube  and  mT 
is  the  total  mass  flow  through  the  entire  inlet  as  defined 
by  equation  3.  This  geometry  is  further  shown  in  Figure 
11-14.  If  the  inlet  lip  splits  a  streamtube  the  values  for 
pressure  recovery  and  mass  flow  are  linearily  interpolated 
to  achieve  the  average  pressure  recovery.  The  values  of 
(pt/ptoo)j  were  taken  from  the  NASA  Ames  printout. 


Outer  Streamline  of  jth  Streamtube 


Capture  streamline 
\ 


Inlet  Lip 


////////  /  //  // 

Body 


Figure  11-14  Geometry  Of  Streamtubes  At  The  Inlet 
Annulus  Used  To  Calculate  Average 
Stagnation  Pressure, 
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III.  CALCULATION  OF  RAMJET  PERFORMANCE 
A.  PERFORMANCE  EQUATIONS 

Ramjet  engine  performance  can  be  quantified  by  the  specif¬ 
ic  thrust,  specific  fuel  consumption,  and  the  coefficient 
of  thrust.  Performance  equations  will  be  developed  for  each 
of  the  above  performance  criterion  which  include  the  concept 
of  additive  drag. 

If  one  dimensional  flow  is  assumed,  the  thrust,  F,  of  a 
ramjet  engine  is  described  by  Netzer  [Ref.  5]  as 


in  v 
e  e 


-  m  v 
o  o 


-  Da  + 


Ae(Pe 


'  V 


(31) 


where  the  subscript  e  refers  to  exit  conditions  and  the  sub¬ 
script  o  refers  to  the  inlet  conditions.  Term  1  is  defined 
as  jet  or  gross  thrust,  term  2  is  the  ram  drag,  term  3  is 
the  additive  drag  and  term  4  is  the  pressure  thrust.  If 
ideal  expansion  is  assumed  (p  =  p  ) ,  equation  1  becomes 

Q  O 


F 


niv 
e  e 


rrtv 
o  o 


D 

a 


(32) 


Now,  defining  the  fuel  air  ratio,  f,  as 


f  »  mf/mo 


(33) 
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the  mass  flow  through  the  ramjet  exit  can  be  written  as 


rii  =  iti  +  =  in  (1  +  f)  (34) 

e  o  f  o 

The  fuel  air  ratio,  f,  is  generally  much  less  than  one. 
Combining  equations  32  and  34  allows  the  total  thrust  to  be 
written  as 


F  =  m  (v 
o  e 


*  v„)  ' 

o  a 


(35) 


Factoring  vq  from  the  first  term  on  the  right  hand  side  gives 


F  '  "ovo  <-T-  -  11 
o 


(36) 


Combining  equation  36  with  the  definition  of  Mach  number  and 
the  speed  of  sound,  aQ,  inserting  an  expression  for  combustor 
energy  balance  and  assuming  that  the  stagnation  pressure 
through  the  inlet,  burner  and  nozzle  are  taken  as  to  be 
constant  allows  thrust  to  be  written  as  a  function  of  fuel 
air  ratio,  heating  value  of  the  fuel  and  the  combustion 
temperature. 


niv 
o  o 


t  r 


l  + 


“ITT 

rr 


-i  i 


p  to 


(37) 
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The  specific  thrust  is  obtained  from  equation  (.37)  by- 
dividing  by  the  mass  flow  rate,  mQ.  Units  of  mQ  are  slug/sec 
and  specific  thrust  is  usually  in  units  of  lb^/ (lbm/sec) . 


SF  =  —  (38) 

mQg 

The  specific  fuel  consumption,  SFC,  is  defined  as  the 
unit  mass  flow  of  fuel  per  hour  per  pound  thrust.  Expressed 
in  the  form  of  an  equation. 

SbOOm^g 

SFC  =  - p-^-  (39) 

In  equation  (39)  m^  has  units  of  slug/sec.  Inserting  equation 
37  yields 


SFC 


3600ni£  g 


A  *  Jsi- 


Vo  1  ’  CT 

p  to 


1]  -  D 


a 


where  m^ 

S 

h 

f 


mass  fuel  flow  rate  (slugs/sec) 
gravitational  constant 
heating  value  of  fuel  (BTU/lb) 
fuel  air  ratio 

specific  heat  capacity  (BTU/lb  °R) 
inlet  stagnation  temperature 


(40) 
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The  thrust  coefficient  is  defined  as 


C 


f 


F 

q  A 

r 


(41) 


where  qQ  is  the  dynamic  pressure  and  A r  is  the  reference 
area.  The  reference  area  is  the  base  area  of  the  five  inch 
projectile.  For  a  ramjet,  mass  flow  at  the  inlet  is  given 
as 


m  =  A  v^o 
o  o  o  o 


(42) 


Combining  equations  (37) ,  (41)  and  (42)  gives 


2A 


ITT 


r2-  [  /  i  +  rr~  -  n  -  rjr 

\  to  qoAr 


(43) 


The  preceding  equations  were  derived  for  a  ramjet  engine 
without  internal  losses.  Losses  in  a  ramjet  engine  that 
should  be  accounted  for  are 

a)  loss  of  stagnation  pressure  in  the  diffuser, 

b)  loss  of  stagnation  pressure  in  the  burner, 

c)  loss  of  stagnation  pressure  in  the  nozzle,  irn 

d)  combustion  efficiency,  nu 


1  ^ 


These  terms  result  in  the  generation  of  a  common  term, 

Sr,  which  does  affect  the  calculation  of  ramjet  performance. 

S  is  defined  as 
r 


‘"b’d’n 


Y-l 

•)T~  -i 


M2 

4  o 


1  + 


nbfh 

CPTto 


1/2 


(44) 


The  square  root  appearing  in  equations  (28),  @0)  and  is  replaced 
by  to  give  ramjet  performance  with  internal  losses  and  with 
an  additive  drag  term.  The  performance  equations  becomes 


SF 


SFC 


["■oVo  (Sr  •  1J  ■  Da1/mos 


3600m£g 

- F - 


2A. 


D. 


csr  -  i)  -  ^ 


(45) 

(46) 

(47) 


o  r 


The  computer  program  developed  by  Fuhs  [Ref.  6]  was 
modified  to  compute  the  additive  drag  term.  Holding  the  values 
of  tt^,  ir^  and  ri^  constant,  the  performance  of  a  ramjet  engine 
with  a  blunt  centerbody  was  predicted.  The  performance  of 
a  ramjet  with  a  conical  spike  inlet  as  predicted  by  the 
program  serves  as  a  baseline  measurement  of  performance. 


IV.  DISCUSSION  OF  RESULTS 


A.  PERFORMANCE  CRITERION 

Performance  of  a  ramjet  engine  can  be  quantified  in 

different  ways.  For  purposes  of  this  report,  the  performance 

was  compared  to  a  ramjet  baseline  configuration  using  a  coni- 

2 

cal  inlet  spike;  the  inlet  capture  area  is  0.0123  ft  [Ref.  7] 
The  performance  of  a  ramjet  with  the  blunt  centerbody 
configuration  was  calculated  for  a  range  of  fuel/air  ratios. 
Performance  is  specified  on  the  basis  of  combusiton  exit 
temperature  (Tt3)  ,  specific  fuel  consumption  and  excess 
thrust.  Of  special  interest  is  performance  when  Cj  equals 
CD.  The  value  of  C^,  which  was  obtained  from  White  [Ref.  8], 
was  0.349. 

The  combustor  exit  temperature  is  calculated  by  the  ram¬ 
jet  performance  program  as  an  output  variable.  Current  tech¬ 
nology  limits  the  steady  combustor  exit  temperature  to  about 
4400°R.  The  value  chosen  for  analysis  purposes  is  4422°R. 
Above  this  temperature,  the  combustor  and  exit  nozzle  will 
melt,  if  run  for  a  continuous  period  of  time.  The  value  of 
Ttj  is  generally  less  than  T^  at  the  stoichiometric  fuel/air 
ratio  for  the  hydrocarbon. 

For  the  final  seconds  of  flight,  it  is  possible  to  boost 
the  fuel/air  ratio  to  the  stoichiometric  value.  The  boost 
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in  fuel/air  ratio  will  increase  the  amount  of  excess  thrust 
available  for  terminal  maneuvering.  Since  the  GLM  will 
explode  at  flight  termination,  the  combustor  exit  temperature 
is  of  little  importance  in  terminal  flight. 

A  further  basis  for  performance  comparison  is  the  effect 

on  relative  detection  range  of  seeker  aperture.  Seeker 

2 

aperture  is  equal  to  irrn  . 

Specific  fuel  consumption  (SFC)  was  chosen  as  a  compari¬ 
son  standard  because  it  provides  a  good  measure  of  ramjet  fuel 
efficiency.  Excess  thrust  is  used  as  a  measure  of  the  ability 
of  the  guided  projectile  to  maneuver  in  the  terminal  phase  of 
the  encounter.  Regardless  of  the  seeker  accuracy,  lack  of 
thrust  above  the  amount  of  flight  drag  will  result  in  the 
inability  to  pursue  a  maneuvering  target.  Excess  thrust 
coefficient,  C£e,  is  calculated  as 


Cfe  CTt3> 


Cf(Tt3)  -  cD 


(48) 


where  is  airframe  drag  and  Cf(Tt3)  is  the  thrust  coeffi¬ 


cient  for  T 


t3 


Two  values  of  are  of  interest.  One  is 


the  maximum  allowed  combustor  exit  temperature  for  steady 
operation  which  is  4422°R.  The  other  is  Tt3  for  stoichiomet¬ 
ric  combustion  which  is  5968°R;  stoichiometric  fuel/air 
ratio  is  0.06. 
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B.  RESULTS  OF  CALCULATIONS 

Figure  IV-1  plots  T^  versus  r^/r^  at  a  value  of  C^  *  C£. 

As  can  be  clearly  seen,  as  the  ratio  of  r./r  increases, 

l  n 

or  as  the  centerbody  shrinks,  T^  decreases.  Tt3  reaches 

the  maximum  sustained  operating  value  of  4422°R  when  r^/rn 

is  1.63;  r,/r  equal  to  1.63  is  the  minimum  value  for  sus- 

tained  operation.  For  stoichiometric  combustion  (f*0.06), 

rT /r  could  possibly  be  reduced  to  a  value  as  small  as  1.50. 
n 

From  the  minimum  value  of  r./r  obtained  in  Figure  IV-1, 

n 

Figure  IV-2  is  entered.  Figure  IV-2  is  a  plot  of  excess 

thrust  coefficient  versus  r, /r  and  of  stoichiometric  excess 

L  n 

thrust  coefficient  versus  r./r  .  Stoichiometric  excess 

l  n 

thrust  is  defined  as  C^  at  stoichiometric  conditions  (f=0.060) 
minus  C^. 

Using  a  value  of  r./r  equal  to  1.63,  the  ramjet  will 
have  excess  maneuvering  thrust  only  at  values  of  r^/rn  greater 
than  1.63.  Use  of  stoichiometric  fuel/air  mixture,  during 
the  last  seconds  of  the  encounter  with  a  maneuvering  target, 
results  in  excess  thrust  being  generated  for  r^/rn  greater 
than  1.51. 

Figure  IV- 3  is  a  graph  of  the  specific  fuel  comsumption 

(SFC)  versus  r./r  .  SFC  decreases  as  r, /r  increases.  The 
l  n  l  n 

baseline  ramjet  has  a  SFC  equal  to  2.02  for  equal  to  C^. 

For  r^/rn  equal  to  1.63,  the  ramjet  with  a  blunt  nose  has 
65%  higher  SFC. 
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Figure  IV-1  Combustor  Exit  Temperature  As  a  Function 
Of  rL/rn  At  Cf  *  CD  »  0.349 


I 


.5 


Cf(5968°R)-CD 


Cf(44Z2»R)-CD 
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Figure  IV-4  is  a  plot  of  relative  detection  range  of  an 
IR  seeker  for  various  rT/r  .  The  signal-to-noise  ratio  of 
the  IR  seeker  and,  therefore,  the  detection  range  is  propor¬ 
tional  to  the  area  of  the  lense.  The  lense  area  is  directly 

related  to  the  square  of  the  lense  radius.  As  r./r  in- 

l  n 
2 

creases,  the  detection  range  decreases  as  ( rn/ rL ' 

Mass  flow  ratio  was  defined  in  the  text  near  Table  II- 3. 

The  ratio  is  interpreted  as  the  fraction  of  the  mass  flow 

actually  captured,  mg,  relative  to  mass  flow  into  an  area 
2 

rcr^.  Figure  IV-5  has  mass  flow  ratio  plotted  as  a  function 

of  r,/r  .  Also  Figure  IV-5  has  a  plot  of  capture  streamtube 

radius  ratio,  r  /r.  ,  as  a  function  of  r./r  .  For  the  mini- 
S  L  l  n 

mum  value  of  r^/rn  equal  to  1.63,  the  capture  radius  ratio 
is  0.84.  Likewise,  the  mass  flow  ratio  is  0.70. 

Figure  IV-6  shows  the  variation  of  additive  drag  coeffi¬ 
cient  with  r./r  .  At  the  minimum  value  of  r./r  of  1.63, 
ij  n  l  n 

Cdad  0*085.  The  airframe  drag  coefficient  is  0.349.  Hence, 

Cdad  *s  24%  of  CB.  To  understand  the  behavior  of  C^a(j  as  a 

function  of  ri/rn  refer  to  Figure  IV- 7.  For  rL/rn  equal  to 

1.1,  which  is  illustrated  in  Figure  IV-7(a),  the  capture 

streamline  between  bow  shock  and  lip  of  the  annulus  is  steep. 

Also,  the  pressure  is  large  near  the  axis.  Consequently, 

Cdad  rs  large  for  small  r^/r^. 

The  greater  the  ratio  of  r^/rn  becomes  the  better  ramjet 

performance  becomes.  However,  to  achieve  even  marginal 

performance  at  r./r_  ■  1.63,  an  unrealistic  penalty  in  the 

l  n 


RELATIVE  DETECTION  RANGE 


Ratio  And  Capture  Radius  Ratio 


Figure  IV-6  Additive  Drag  Coefficient  As  A  Function  Of 
rL/rn  For  An  Inlet  At  The  Shoulder  (Z=1.00) 
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1.4 


b)  rL/rn  * 

Figure  IV- 7  Bow  Shock  Wave  And  Capture  Streamtube 
For  r./r  Equal  To  1.1  And  1.4 


I  -TF 


form  of  a  63%  loss  of  relative  detection  range  must  be  paid. 
Lower  values  of  r./r_  are  unable  to  be  obtained  in  a  thrust 
equals  drag  configuration  due  to  high  values  of  additive 

I 

drag  and  poor  ir^. 

Table  IV-1  contrasts  ramjets  with  the  spike  inlet  and 
two  blunt  nose  inlets.  All  three  ramjets  have  identical 
thrust  coefficients  equal  to  airframe  drag  coefficient.  Due 
to  large  specific  thrust,  the  ramjet  with  spike  inlet  re¬ 
quires  considerably  less  mass  flow  rate.  The  ramjet  with 
spike  inlet  has  superior  specific  fuel  consumption.  Excess 
thrust  is  comparable  for  the  three  ramjets. 
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Table  VI-1 

Comparison  Of  Spike  Inlet  Versus  Two  Blunt  Nose  Inlets 


Ramjet  Using 

Parameter 

Spike 

Inlet 

Blunt  Nose  Inlet 

Best  Ramjet  * 

Best  IR  * 

Airframe  drag 
coefficient,  CD 

0.349 

0.349 

0.349 

Additive  drag 

coefficient,  C^a^ 

0 

0.018 

0.085 

Thrust  coefficient, 

0.349 

0.349 

0.349 

Combustor  exit  temperature 
°R,  Tt3  4122 

3284 

4422 

Mass  flow  rate,  lb  /sec 

m 

9.5 

18.3 

13.6 

Specific  Thrust, 
lbf/(lbm/sec) 

67.2 

34.8 

46.8 

Specific  fuel  consumption 

Clbm/hr) / 

*2.02 

2.68 

3.35 

Maximum  thrust  coefficient 

Cf,  at  stoichiometric 

fuel  air  ratio  0.505 

0.798 

0.509 

Combustor  exit  tempera¬ 
ture,  Tt3,  °R 

5698 

5698 

5698 

Relative  detection 
range 

0.33 

0.29 

0.37 

*  Best  ramjet  is  obtained  for  r./r  equal  to  1.854  at  least 

for  the  range  of  r./r  investigaied  here.  Best  IR 

detection  capability  occurs  for  rT /r  equal  1.63. 

l  n 
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V.  CONCLUSIONS 


Conversion  of  the  Navy  SALGP  from  a  rocket  to  a  ramjet 
severely  degrades  detection  capability,  if  an  axisymmetric 
inlet  is  used.  Note  the  values  of  0.33,  0.29,  and  0.37  for 
relative  detection  range  from  Table  IV- 1. 

The  ramjets  with  blunt  nose  inlets  suffer  in  both  speci¬ 
fic  thrust  and  specific  fuel  consumption.  The  poor  perfor¬ 
mance  is  due  to  two  causes.  First,  the  high  value  of  ^ 
and  second,  poor  pressure  recovery  of  the  inlet  leading  to 
poor  performance. 

Surprisingly,  the  excess  thrust  coefficient  of  ramjets 
with  blunt  noses  is  very  competitive  with  the  ramjet  with 
spike  inlet. 

In  view  of  the  foregoing  conclusions,  pop-out  or  retrac¬ 
table  scoop  inlets  appear  much  more  attractive  in  spite  of 
structural  constraints  and  space  considerations.  Use  of  an 
axisymmetric  inlet  located  in  the  nose  of  a  5-inch  projectile 
is  a  severe  design  compromise. 
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APPENDIX  A 


NASA  AMES  COMPUTER  PROGRAM  DESCRIPTION 

I.  INTRODUCTION 

Two  different  computer  programs  were  used  to  model  the 
flow  arouiii  the  blunt  nosed  body.  The  coordinate  system 
used  in  the  thesis  and  the  NASA  Ames  computer  program  is 
shown  in  Figure  I I -6. 

The  first,  IMPLCBO,  is  a  modified  version  of  AXI -BLUNT; 
see  Kutler,  Chakravarthy ,  and  Lombard  [Ref.  9] .  Programmed 
to  predict  the  supersonic  flow  over  a  three  dimensional 
body,  the  algorithm  predicts  shock  shape  and  location  as  well 
as  flow  parameters  at  equally  spaced  interior  points  between 
the  body  and  the  shock. 

The  second  program,  OGIVE,  solves  for  the  supersonic  flow 
parameters  around  a  three  dimensional  wing  body  configuration; 
refer  to  Kutler,  Reinhardt,  and  Warming  [Ref.  10].  Originally 
designed  for  flow  prediction  around  a  delta  winged  spacecraft, 
the  program  will  provide  flow  field  parameters  at  various  radii 
from  the  centerbody  as  well  as  at  various  locations  along 
the  body. 

II.  PROGRAM  DESCRIPTION 

The  IMPLCBO  program  uses  an  unsteady,  implicit  numerical 
procedure  to  determine  the  supersonic  flow  around  the  body. 
Further  program  description  is  contained  in  Reference  9. 


* 
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The  program  output  is  referenced  to  various  free  stream 
values.  These  values  are  shown  in  Table  A-l. 

The  program  utilizes  the  body  configuration  and  predicted 
shock  points  as  boundary  values  and  solves  the  inviscid 
Euler  or  Navier-Stokes  equations  for  the  interior  points, 
after  a  coordinate  transform. 

The  different  angular  relationships  at  the  shock  front 
are  shown  in  Figure  A-l. 

The  OGIVE  program  utilizes  the  flow  field  parameters 
determined  by  IMPLCBO,  at  9  *  90°  and  X  *  1.00,  as  an  initial 
starting  point  and  determines  flow  field  values  along  the 
cylinder  body  by  a  finite  difference  algorithm.  The  program 
is  further  described  in  Reference  10.  The  program  reference 
values  are  tabulated  in  Table  A-l. 

The  program  describes  the  body  shapes  by  analytical  approx¬ 
imations.  These  approximations  consist  of  coordinate  posi¬ 
tions  and  the  slope  of  the  body  contour  as  a  function  of  the 
distance  along  the  body  axis. 


Table  A-l 

Reference  Values  For  OGIVE  And  IMPLCBO  Programs 


IMPLCBO  (region  1) 


output  variable 


reference  value 


P/PINF 

RHO/RINF 

U/QINF 

V/QINF 


Poo 


/y  M^;  at  >1^*3. 0 
/yM  ;  at  M  =3.0 

OO  QO 


OGIVE  (region  2) 


R 

r  -  1.0 

T1 

P 

Pt°°  ; 

P»/Ptco*0-0272237  at 
M^-3.0 

RHO 

Pt-  ; 

p  /p..  =0.0762263  at 
M.-3.0 

QINF 

<lco“Vco/V; 

’q«»[l  +  5/M f] 1/2 

U 

^00 

Z  component  of  flow 
velocity 

V 

R  component  of  flow 
velocity 

W 

0  component  of  flow 
velocity 

Z 

rn  ’ 

length  along 
cylinder  body 

Figure  A-l  Angular  Relationship  At  The 
Shock  Front  For  5,  S  ,  0 
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APPENDIX  B 


MASS  CONTINUITY  METHOD  OF  DETERMINING  STREAMLINE  SHAPE 

A  computer  program  has  been  written  for  the  HP  9830  to 
calculate  the  capture  streamline.  The  program  is  designed 
to  reduce  the  output  data  obtained  from  the  NASA  Ames  com¬ 
puter  output  and  compute  the  streamtube  position  at  various 
points  in  the  flow  field  around  the  blunt  nosed  body. 

The  streamtube  position  is  determined  by  comparing  the 

mass  flow  of  air  entering  the  bow  shock  to  the  integral  of 

mass  flux.  The  streamtube  radius  at  the  bow  shock,  r  ,  is 

s 

determined  by  comparing  the  mass  flow  at  the  inlet  to  mass 
flow  at  the  nose. 

The  value  of  the  integral  in  equation  (22)  or  (23)  is 
determined  by  assuming  a  linear  change  in  pressure,  density 
and  velocity  between  successive  data  points.  The  program 
interpolates  to  find  the  distance  from  the  axis  that  the 
streamtube  must  be  to  satisfy  the  mass  flow  requirements  of 
the  inlet.  Values  of  X  and  Y  as  shown  in  Figure  I I- 6  are 
calculated,  and  an  interpolated  pressure  is  determined. 

The  value  of  the  integral  in  equation  (22)  and  (23)  is  cal¬ 
culated  using  the  rectangular  method.  Equations  (22)  and 
(23)  are  developed  in  Chapter  III  of  this  thesis. 
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Table  B-l 

Variables  For  Capture  Streamtube  Computer  Program 


Symbol  for 
HP  9830 

AO 

A7 

C 

Cl 

C2 

F 

G 

1 

P7 

S7 

W 

Q9 

R9 

S 

X7 

Y7 


Definition 

value  of  integral  in  equation  (22)  or  (23) 
between  point  i  and  i+1 

remaining  area  between  tabulated  values 

squared  radius  of  streamtube  at  shock  front 
(square  inches) 

constant  for  integral  when  in  region  I 

constant  for  integral  when  in  region  II 

percentage  of  distance  Y7  is  between  Y(I) 
and  Y(I+1) 

dummy  variable  for  area  calculation 
indexing  variable 

interpolated  value  of  pressure  X7,Y7 
slope  of  area  curve 

dummy  variable  used  in  calculation  of  C 
value  of  Q^,  equal  to  (1+  (5/M^))'1^2 
value  of  oj ptoo  at 

summation  of  integral  area  to  the  ith  point 
(inches  squared) 

interpolated  value  of  X,  determined  from  Y7 
(inches) 

interpolated  value  of  Y  (inches) 


Input  variables 

A  angular  value  of  data  points  above  projectile 

centerline  (degrees) 

19  number  of  data  points  in  data  file 

Z  location  of  data  points  in  relation  to  the 

projectile  shoulder.  1*  in  front  of  shoulder 
2*  on  the  body  cylinder 
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Table  B-l  con't 


Symbol  for  Definition 

HP  9830 

Dimensioned  Variables 

G(I)  integral  value  between  i  and  i+1 

L(I)  integral  value  at  ith  point 

P(I)  pressure  ratio  at  ith  point 

R(I)  density  ratio  at  ith  point 

U(I)  component  of  flow  velocity  along  the  body 

in  the  Z-direction 

V(I)  component  of  flow  velocity  along  the  body 

in  the  r  direction 

Y(I)  distance  from  the  body  centerline  to  the 

data  point 


Table  B-2  Capture  Streamtube  Program  Listing 
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APPENDIX  C 


ADDITIVE  DRAG  COMPUTATION 

The  program  is  designed  to  calculate  the  additive  drag 
and  coefficient  of  additive  drag  for  an  inlet.  The  reference 
area  is  the  projectile  base  area.  The  program  is  written  in 
BASIC  for  the  HP  9830  desk  top  calculator. 

The  input  data  consists  of  the  point  position  (X  and  Y) 
and  the  static  pressure  at  the  point.  The  three  variables 
X,  Y,  and  p  are  determined  by  first  running  the  program 
described  in  Appendix  B  of  this  thesis. 

Equations  (27)  through  (30)  used  in  this  calculation  are 
developed  in  Section  III  of  this  thesis. 

The  computer  symbols  are  listed  and  defined  in  Table  C-l 
and  the  program  listing  is  in  Table  C-2. 


Table  C-l 

Variables  For  Coefficient  Of  Additive  Drag  Computer  Program 


Variables  Definition 


Reference  area  (square  inches) 

Coefficient  of  additive  drag  between  I  and  1+1 
data  point  (dimensionless) 

total  coefficient  of  additive  drag  (dimensionless) 
total  drag  of  body 

distance  along  body,  Z,  between  data  points 
(inches) 

counting  variable 
counting  variable 

Ratio  of  circle  circumference  to  diameter 
p  /p..  at  input  M 

Distance  from  centerline,  Y  (inches) 

Slant  distance  between  data  points  (inches) 

Surface  area  of  cone  between  data  points 
(square  inches) 

Incremental  drag  between  data  points  (lb) 


Input  variables 


Mach  number  (dimensionless) 
Number  of  data  points 


Dimensioned  variables 


A(I)  Angle  of  streamtube  between  data  points  (degrees) 

P(I)  Pressure  at  Ith  data  point  (dimensionless) 

R(I)  Average  pressure  between  data  points 

(dimensionless) 

X(I)  X  position  of  Ith  data  point  (inches) 

Y(I)  Y  position  of  Ith  data  point  (inches) 


Table  C-2  Coefficient  Of  Additive  Drag  Program  Listing 
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